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ABSTRACT 



We present new thermal IR photometry and spectral energy distributions (SEDs) of eight classical Cepheids (type I) and three type II 
Cepheids, using VISIR thermal IR photometric measurements, supplemented with literature data. We used the BURST mode of 
the instrument to get diffraction-limited images at 8.59, 11.25 and I1.85yum. The SEDs show a IR excess at wavelengths longer 
than lOyum in ten of the eleven stars. We tentatively attribute these excesses to circumstellar emission created by mass loss from 
the Cepheids. With some hypotheses for the dust composition, we estimated a total mass of the envelope ranging from lO""' to 
10"* Mq. We also detect a spatially extended emission around AX Cir, X Sgr, W Sgr, Y Oph and U Car while we do not resolve the 
circumstellar envelope (CSE) for the other stars. The averaged circumstellar envelope brightnesses relative to the stellar photosphere 
areo-CAXCir) = 13.8 + 2.5 %, o-CX Sgr) = 7.9± 1.4 %, o-CW Sgr) = 3.8 + 0.6%,q'(Y Oph) = 15.1 ± 1.4 % and o-CU Car) = 16.3±1.4% 
at 8.59//m. With this study, we extend the number of classical Cepheids with detected CSEs from 9 to 14, confirming that at least a 
large fraction of all Cepheids are experiencing significant mass loss. The presence of these CSEs may also impact the future use of 
Cepheids as standard candles at near and thermal infrared wavelengths. 

Key words, stars: circumstellar matter - stars: variables: Cepheids - stars: mass-loss - stars: imaging - Infrared: stars 



1. Introduction 



IKervella et al.l ( |2006|) discovered in the and K band a circum- 
stellar envelope (CSE) around { Car using the MIDI and VINCI 
instruments from the VLTI. Its typical size is 2 - 3 with a 
contribution of 4 % to the total flux in K. Similar i nterferometric 
detec t ions w ere then reported for other Cepheids (iMerand et alJ 
l2006ll2007h leading to the conclusion that a significant fraction 
of all Cepheids are surrounded by a CSE . From mid- and far- 
IR observations with the Spitzer telescope. iBarmbv et al.l (1201 ih 
also detected extended emission around a significant fraction of 
their sample of 29 classical Cepheids. The case of 6 Cep and it s 
extended emission is well discussed in lMarengo et al.l (l2010bl) . 
These CSEs have an effect on the infrared surface brightness 
technique (IRSB) since the Cepheid appears brighter and also 
on interferometric measurements since the star appears larger 
than it really is. It is therefore necessary to quantify this excess 
(linked to an extended emission) in order to estimate the bias on 
the distance. This is particularly important for the James Webb 
Space Telescope (JWST) that will be able to observe Cepheids 
in distant galaxies. Their distances could be estimated via the 
use of IR P-L relations but the presence of infrared excess will 
degrade the distance accuracy. 

This circumstellar material is also important in the context 
of Cepheid mass-loss since it may play a signi ficant role in 
the pr oblem of the Cepheid mass discrepancy dNeilson et alj 
1201 OaL ...). The infrared excess could be linked to a mass-loss 
mechanism generated by shocks between different layers of the 
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Cepheid's atmosphere during the pulsation cycle. A correla- 
tion between the perio d and the envelope fl ux (relatively to the 
star) was proposed bv iMerand et al.l (|2007|) . A Cepheid with a 
larger pulsation period would have a larger IR excess. On the 
other hand, from pho tometry on a larger sample of Cepheids, 
iNeilson et al] ( 1201 Obi) reach a different conclusion. It is thus es- 
sential to study these CSEs to quantify their contribution and to 
understand how they form. 

To progress on this question, we present new observations 
of a few classical (type I) and type II Cepheids with the VISIR 
instrument of the VLT. The mid-IR wavelength coverage of this 
instrument is well suited to study the infrared contribution of 
the CSEs. The instrument configurations and the data reduction 
methods are detailed in Sect. |2l In Sect. [3] we present the data 
analysis using aperture photometry applied to our images and 
we study the spectral energy distribution of our sample of stars. 
We also look for a spatially resolved component using a Fourier 
technique analysis. Finally we discuss our results in Sect.|4] 

2. Observations and data processing 

The selected sample of classical Cepheids was chosen for their 
brightness, the range of pulsation period (from 4 to 45 days) and 
their angular size so that they are resolvable by long-baseline 
interferometry. Y Oph s eems to exhibit the brightest CSE de- 
tected around a Cepheid dMerand et alJl2007l) and the goal is to 
explore larger distances from the Cepheid. The type II Cepheids 
we selected have been extensively studied and are known to have 
strong IR excesses linked to their CSEs. They will be used in our 
analysis both to compare the properties of their CSEs with those 
of classical Cepheids, and to validate our modeling approach to 
estimate the thermal infrared excess (based on the results of pre- 
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vious works on type II Cepheids). We present in Table [T] some 
relevant parameters for our sample of observed Cepheids. 

The observations were performed using the ESO mid- 
infrared instrument VISIR installed on UT3 at the VLT (Paranal, 
Chile). This instrument works in the (8 - 13//m) and Q 
(16 - 24jum) atmospheric windows and provides imaging and 
long-slit spectroscopy. As we need the best sp atial resolution, we 
chose the BURST mode (iDoucet et al.ll2006l) of VISIR to over- 
come the atmospheric seeing. In this mode, short exposure im- 
ages (i.e. faster than the atmospheric coherence time) are taken 
in order to freeze the turbulence. Each short exposure frames so 
presents one principal speckle. With sets of thousands of short 
exposure frames we can apply a shift-and-add process. This en- 
hances the quality of the data and enabled us to reach the diffrac- 
tion limit of the telescope. To correct for instrumental artefacts 
and background thermal emission, we applied the classical chop- 
nod technique. 

The raw data of the VISIR BURST mode are data cubes in 
which the star appears in a different position in the field (the 
chop-nod positions, perpendicular each other) for each frames. 
Our first step consisted of a classical subtraction of the chopped 
and nodded images in order to remove the thermal background 
and in storing them in data cubes that contain thousands of 
frames (~ 20000). To have the best diffraction-limited images, 
we selected 50 % of the best frames according to the brightest 
pixel (as tracer of the Strehl ratio). We then proceeded to a pre- 
centering (at a integer pixel level), a spatial resampling by a fac- 
tor of 4 using a cubic spline interpolation and a fine recentering 
using a Gaussian fitting (at a precision level of a few milliarcsec- 
onds). The resulting cubes were then averaged to get the final im- 
age used in the data analysis process (see Sect.O. This raw data 
processin g has already been used and has proven its efficiency 
(see e.g. 'Kervella & Domi ciano de Souzall2007t iKervella et alJ 
l2009l: |Gallenne et al. 2011). 

The observations were carried out on the nights of 2008 
May 23-24. Table |2] lists the sequence of our observations in- 
cluding the reference stars that were observed immediately be- 
fore and after the Cepheids in the same instrumental setup. 
Series of observations were obtained in three filters; PAHl, 
PAH2 and SiC (respectively 8.59 + 0.42 /zm, 11.25 + 0.59 |/m 
and 11.85 + 2.34 yum). We chose the smallest pixel scale of 
75 mas/pixel (before resampling) to ensure a proper sampling of 
the FWHM of the Airy pattern (average FWHM ~ 4 pixels for 
the reference stars). For an unknown reason, the nodding posi- 
tion was out of the detector for the observations #15, #17, #19 to 
#22, #25 to #28 and #31 to #34 but we still have enough frames 
(~ 10000) to correctly subtract the background emission. The 
observations #23 and #24 were chopped and nodded out of the 
detector and were not included in the analysis. Due to the lower 
sensitivity of VISIR in the SiC filter, we could not recenter the 
individual VISIR images for the observations #40 and #44. The 
#4 sequence was not included because of very low signal, maybe 
caused by some clouds, leading to a bad recentering of the indi- 
vidual images. 
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Fig. 1. Relative atmospheric transmission for our observing 
nights. Both nights are normalised to the mean value over all 
calibrators observed. The dotted lines denote a 5 % variation. 



a Fourier technique analysis on our average VISIR images we 
also search for the presence of extended emission. 

3.1. Atmospheric evolution 

We studied the sky transparency evolution using o ur photomet- 
ric ref erence stars. The photometric templates from lCohen et al.l 
([T999I) were used to plot the temporal relative atmospheric trans- 
mission (relative to the average value of all calibrators observed 
for each night). We see in Fig. [T]that the atmosphere was pretty 
stable during the observing nights. In the first night, the relative 
standard deviation is below 5 % while it is only 3 % in the second 
night. 

3.2. Pliotometry 

We carried out a classical aperture photometry to assess the flux 
density for each st ar in each filter Photometric templates from 
ICohenetalJ 119991) were used to have an absolute calibration of 
the flux density taking into account the filter transmissiorQ using: 



F(ref) = 



where is the reference irradiance (from the Cohen et al. tem- 
plates), Ta the filter transmission and AA the filter bandpass. 
We applied an airmass correction factor; 

/^coiT = /^obsXC(/l,AM) 

with C(A,AM) taken from lSchutz & SterzikI (|2005|) : 



3. Data analysis 

In this section we first study the evolution of the atmospheric 
properties (seeing, transparency) during the observations for 
each night. We then present the aperture photometry we ap- 
plied to our average images. We then combine these measure- 
ments with other photometric data we retrieve from the litera- 
ture to study the spectral energy distribution of each star. Using 



C(A,AM) = 1 -I- 



0.104 

0.220 — (^-8.6//m) 



(AM - 1) 



No aperture correction is required since we kept the same, 
relatively large, aperture radius (1.3") for both the Cepheids and 



' Filter transmission profiles are available on the ESO website 
http://www.eso.org/sci/facilities/paranal/instruments/visir/inst/index.htmll 
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Table 1. Some relevant parameters of our Cepheids. 



Stars 


pa 


MJDo" 


< V >^ 


<K>'' 








Type 




(days) 








(mas) 


(mas) 


FFAql 


4.4709 


41575.928 


5.37 


3.49 


0.88 + 


0.05 


2.81 ± 


0.18 I 


AXCir 


5.2733 


51646.100 


5.88 


3.76 


0.70 ± 


0.06 


3.22 ± 


1.22 I 


XSgr 


7.0128 


51653.060 


4.55 


2.56 


1.47 + 


0.04 


3.00 ± 


0.18 I 


Aql 


lAlbl 


36084.156 


3.90 


1.98 


1.84 ± 


0.03 


2.78 ± 


0.91 I 


WSgr 


7.5950 


51652.610 


4.67 


2.80 


1.31 ± 


0.03 


2.28 ± 


0.20 I 


YOph 


17.1242 


51652.820 


6.17 


2.69 


1.44 ± 


0.04 


2.04 ± 


0.08 I 


UCar 


38.8124 


51639.740 


6.29 


3.52 


0.94 + 


0.06 


2.01 ± 


0.40 I 


SV Vul 


45.0121 


43086.390 


7.22 


3.93 


0.80 + 


0.05 


0.79 ± 


0.74 I 


RSct 


146.50 


44871.500 


6.70 


2.27 






2.32 ± 


0.82 II 


AC Her 


75.010 


35097.300 


7.90 


5.01 






0.70 ± 


1.09 II 


K Pav 


9.0814 


46683.569 


4.35 


2.79 


1.17 + 


0.05 


6.00 ± 


0.67 II 



Notes. P is the pulsation period and MJDq denotes the reference epoch (MJDo = JDo - 2400000.5). <V > and < K > are the mean apparent V and 
K magnitudes, the angular diameter and n the trigonometric parallax. Type denotes the classical (type I) or type II Cepheids. 



Feast et al.l d2008h for k Pav ; from lSamus et alj ( l2009h for the others. 

Samus et al J j2009l) for SV Vul, FF Aql, rj Aql, R Set and AC Her ; from lFeast etal] ( [20081) for K Pav ; from lBerdnikov & CaldwelH 



from 

from 

( 1200 Ih for the others.^ 

f"^' fromlpernie et al. (1 19951) for the classical Cepheids ; from lSamus et al.l i2009l) for the type II. 

'■^^ fromlw elch et al. ( 1984) for FF Aql and X Sgr ; from DE NIS for AX Cir andW Sgr ; from Barne s et alj (|199 7|) for ?/ Aql ; frorn 

iLanev &"stoble ( 1992) for Y Oph, U Car and S V Vul ; froml Taranova et al.l ( 120101) for R Set and AC Her ; fromlpeast et alj (120081) for k Pav. 

limb-darkened angula r diameters from[ Kervella et al. ( 20041^) for X Sgr, ?? Aql, W Sgr a nd Y Oph ; fro mlGroenewegenT i 2007|) for FF Aql ; 
predict ed diameter fr om'Mosk alik & Gorvnva C2005i) for AX Cir ; fromlFea st et al.l (l 2008h fo r k Pav ; fromlGr p eneweg M (l2008l) fo r the others. 

fr omlBenedict et a l. ( 2007 ) for FF Aql, X Sgr and W Sgr ; from lHoffleit & Jaschekl j 19911) for U Car ; from lMerand et alj ( 120071) for Y Oph; 
from lPerrvman et al.l (il997i) for the others. 



their respective reference stars. We then estimated the flux den- 
sity of our Cepheid samples; 

7^/ \ ^corrCcep) 
^^corr(ref) 

where cep and ref stand for the Cepheid and reference stars. 

The measured flux densities are summarized in Table |3] The 
uncertainties include the statistical dispersion of the aperture 
photometry, the dispersion of the calibrator flux densities over 
the night, and the absolute calibration uncertainty. Some final 
averaged images of our sample are presented in Fig. [14] The vis- 
ibility of several diffraction rings reflects the quality of the data. 

3.3. Spectral energy distribution of classical Cepheids 

We have collected additional photometric measurements in the 
literature from 0.4 yum to 100 //m to analyse the spectral energy 
distribution of our Cepheids. As they are pulsating stars, the 
SEDs vary during the pulsation cycle and we have to take into 
account the phase of pulsation when retrieving the data. 

To estimate the magnitudes at our phase of pulsation we re- 
trieved light curves from the literature when available, that we 
plotted as a function of phase (all phases were computed with 
the ephemerides from Table [Tjl. We then applied th e Fourier de- 
composition technique (see e.g. lNgeow et al.|[2003h ; 



= flo + ^ fl/ cos(2m(fi + bi) 



where a, and bj are the Fourier amplitudes and phases for the or- 
der i to fit and (f> is the pulsation phase. The parameter n depends 
on the number of data points and on the light curve amplitude 
(higher amplitude require higher order fit). We used n - 4 for 
all stars. Examples are shown in Fig.[2]for two stars. The coeffi- 
cients of the fit were then used to compute values at our phases. 



We used this method for the Cepheids FF Aql (in B, V and R 
bands), AX Cir {B, V), X Sgr (B, V, J, H, K), jj Aql (B, V), W Sgr 
(B, V), Y Oph (B, V), U Car (B, V) and SV Vul (B, V). When 
data points were not equally spaced in phase, the direct Fourier 
fit produced unrealistic oscillations, we so used a periodic cu- 
bic spline interpolation. To evaluate the uncertainty we used the 
total standard deviation of the residual values (bottom panels of 
Fig. [3. The uncertainty is considered the same for all phases. 
We used this kind of interpolation for the stars rj Aql (in J, H, K 
bands), Y Oph {J, H, K), U Car (7, H, K), S V Vul (J, H, K) and 
Pav {V,J,H,K). 

However there is a lack of light curves in the literature for 
some wavelengths and it is not possible to correct all the data for 
the phase mismatch. We then chose to use the amplitude of the 
light curves A,i as additional uncertainties on the magnitude due 
to the phase mismatch. In addition we know that A,i is decreasing 
with w avelength. Using 5 1 Galactic Cepheids, iLanev & Stobid 
(1 19931) plotted the amplitude of the J,H and K light curves 
as a function of the pulsation period (their Fig. 8). There are 
mainly two regions; 0.5 < logP < 1.0 with Aj ~ 0.1 mag 
and log P > 1 .0 with Aj ~ 0.2 mag and these values decrease 
with the wavelength {Ak ~ 0.08 mag for logP < 1.0 and 
Ak ~ 0.15 mag for logP > 1.0). Therefore when the collected 
additional data do not correspond to same phases, we chose as 
additional uncertainties on the magnitude due to the phase mis- 
match the values 0.1 mag if 0.5 < logP < 1.0 or 0.2 mag if 
logP > 1.0 for 1 < A < 3.5 yum and for A > 3.5 ;um we chose 
0.05 mag for all periods. We are conscious that this can over- 
estimate the uncertainties but this will enable us to completely 
separate an IR excess detection from a magnitude error on the 
SED curves. On the other hand this can also prevent the detec- 
tion of a small IR excess. 

Some of our targets were observed on two different nights 
and they can have different phases (especially the short periods). 
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Table 3. Measured flux densities of our Cepheid sample. 



Stars 


MJD 


Filter 


Flux density 




Flux density 


Excess 








(W/mV/um 






(Jy) 


(%) 


FF Aql 


54 611.403 


PAHl 


9.23 


+ 


0.25 X 


10" 


14 


2.27 ± 0.06 


1.6 : 


t 2.8 




54611.411 


SiC 


2.79 


+ 


0.08 X 


10" 


4 


1.28 ± 0.06 


3.7 : 


t 3.0 


AXCir 


54 611.076 


PAHl 


6.97 


+ 


0.19 X 


10" 


14 


1.72 ± 0.05 


-0.8 


+ 2.9 


XSgr 


54 610.104 


PAHl 


2.31 


+ 


0.11 X 


10- 


li 


5.69 ± 0.27 


14.1 


± 5.4 




54611.112 


PAHl 


2.11 


+ 


0.10 X 


10- 


3 


5.20 ± 0.25 


3.8 : 


t 9.8 




54610.111 


PAH2 


7.82 


+ 


0.20 X 


10- 


14 


3.31 ± 0.08 


11.9 


± 2.9 




54611.119 


SiC 


6.27 


+ 


0.15 X 


10- 


4 


2.89 ± 0.07 


5.6 : 


t 8.1 


77 Aql 


54 610.282 


PAHl 


3.96 


± 


0.14 X 


10- 


i 


9.73 ± 0.35 


0.38 


± 3.6 




54 611.248 


PAHl 


4.30 


+ 


0.16 X 


10- 


3 


10.6 ± 0.39 


9.0 : 


t 4.1 




54 610.289 


PAH2 


1.39 


+ 


0.05 X 


10- 


3 


5.89 ± 0.21 


2.6 : 


t 3.7 




54 611.255 


SiC 


1.28 


+ 


0.05 X 


10- 


13 


5.90 ± 0.23 


9.6 : 


t 4.3 


WSgr 


54 610.213 


PAHl 


1.70 


+ 


0.07 X 


10- 


— 


4.19 ± 0.17 


10.4 


± 4.6 




54 611.131 


PAHl 


1.82 


+ 


0.05 X 


10- 


13 


4.48 ± 0.12 


18.2 


± 3.3 




54 610.220 


PAH2 


6.16 


+ 


0.36 X 


10- 


4 


2.61 ± 0.15 


16.5 


± 6.8 




54611.139 


SiC 


5.30 


+ 


0.13 X 


10- 


13 


24.4 + 0.6 


16.2 


± 2.9 


YOph 


54 610.126 


PAHl 


2.01 


+ 


0.05 X 


10- 


13 


4.95 ± 0.12 


6.8 : 


t 2.7 




54 611.170 


PAHl 


2.00 


+ 


0.05 X 


10- 


13 


4.93 ± 0.12 


6.3 : 


t 2.7 




54 610.134 


PAH2 


6.62 


+ 


0.25 X 


10- 


14 


2.80 ±0.11 


2.4 : 


t 3.9 




54 611.177 


SiC 


5.83 


+ 


0.15 X 


10- 


14 


2.69 ± 0.07 


4.6 : 


t 2.7 


UCar 


54 610.546 


PAHl 


1.07 


+ 


0.02 X 


10- 


13 


2.64 ± 0.05 


32.1 


± 2.5 




54 611.035 


PAHl 


1.06 


+ 


0.05 X 


10- 


14 


2.61 ± 0.12 


30.9 


± 6.2 




54 611.042 


SiC 


2.91 


± 


0.2 X lO^'* 


1.34 ± 0.09 


20.9 


± 8.3 


SV Vul 


54 611.365 


PAHl 


8.62 


+ 


0.17 X 


10- 


14 


2.12 ± 0.04 


25.0 


± 2.5 




54 611.372 


SiC 


2.35 


+ 


0.07 X 


10- 


14 


1.08 ± 0.03 


15.1 


± 3.4 


RSct 


54 610.236 


PAHl 


6.74 


+ 


0.16 X 


10- 


3 


16.6 ± 0.4 


134 


± 6 




54 611.189 


PAHl 


7.05 


+ 


0.17 X 


10- 


13 


17.4 ± 0.4 


145 


± 6 




54 610.243 


PAH2 


2.56 


+ 


0.06 X 


10- 


13 


10.8 ± 0.3 


153 


± 7 




54 611.196 


SiC 


2.24 


+ 


0.06 X 


10- 


13 


10.3 + 0.3 


158 


+ 7 


AC Her 


54 610.349 


PAHl 


9.18 


+ 


0.25 X 


10- 


13 


22.6 ± 0.6 


8273 


± 228 




54 611.340 


PAHl 


9.40 


+ 


0.23 X 


10- 


13 


23.2 ± 0.6 


8473 


± 210 




54 610.357 


PAH2 


10.1 


+ 


0.28 X 


10- 


13 


42.8 ± 1.2 


26705 


± 743 




54 611.347 


SiC 


8.33 


+ 


0.21 X 


10- 


13 


38.4 ± 1.0 


25586 


+ 647 




54 611.227 


PAHl 


1.93 


+ 


0.05 X 


10- 


13 


4.75 ± 0.12 


22.2 


± 3.2 




54 611.235 


SiC 


5.39 


+ 


0.14 X 


10- 


14 


2.48 ± 0.06 


15.3 


± 3.0 



Notes. The flux density was measured over an aperture radius of 1.3". The parameter a is the measured flux density excess. 



We so considered the SED according to an intermediate phase 
between the observations. 

The photospheric emission was modelled with tabulated stel- 
lar atmosphere models obtained with the ATLAS9 simulation 
code from ICastelli & Kurucz (2003). We have chosen a grid 
which was computed for solar metallicity and a turbulence ve- 
locity of 2 km s-'. We then interpolated this grid in order to com- 
pute spectra for any effective temperature and any surface grav- 
ity. The spectrum was multiplied by the solid angle of the stellar 
photosphere, 7t0^^/4, where 0ld is the limb-darkened angular 
diameter. We adjusted the photometric data to the model tak- 
ing into account the spectral response of each instrument. We 
assume that there is no detectable excess (5 % or less) below 
2.2 jum and all the photometric measurements bluer than the K 
band are used to fit the angular diameter and the effective temper- 
ature. However we should note that Teff and 0]^u are correlated 
variables in this fit. We did not adjust the surface gravity, since 
the broadband photometry is mostly insensitive to this parameter 
and its value was also chosen in the Uterature. 

All flux densities < 3 fim are corrected for interstellar extinc- 
tion A,i = R, ] E{B - V) us i ng the total-to-selective absorption 
ratios R,i from lFouque et alj ( l2003h and Hindsley & Befl ( 198§, 
and the color excess E(B - V) from lFouque et alJ (I2007i) . Fluxes 



in any other longer wavelengths are not corrected for the inter- 
stellar extinction, which we assume to be negligible. 



3.3.1. FF Aql 



We selected log^ - 2.05 fromlLucketaD (l2008l hereafter|LOl 
as fixed parameter. The adjusted SED is presented in Fig. [3] The 



V a nd R valu es we re ta ken from light curves of 



Berdnikovl 



(120081 hereafter IbOSI) and iMoffett & BarnesI (11984 ). hereafter 
|M84 . The J, H and photometrv are from lWelch e t al. ( 198|) 



and correspond to the mean values (there was not a good cov- 
erage in phase to estimate accurate light curves, we so de- 
cided to take the mean values and their standard deviations). We 
also added photometric values from the Infrared Array Camera 
(IRAC: 3.6, 4.5, 5.8 and 8/im) and the Multiband Imaging 
Photo meter {MIPS: 24 //m) in stalled in t he Spi tzer space tele- 
scope dMarengo et al.l l20I0aL hereafter IMIOI) . We also use 
broadband photometr y from the Infrared Astr onomical Satellite 
{IRAS: 12 and 25 um. iHelou & Walkeil[T988h . fr om flie AKARI 
satellite IRC point source catalogue (9 and 1 8 fim. Ilshihara et aP 
1201 Oh and from the Mi dcourse Space Experiment {MSX: 8.2 8, 
12.13, 14.65. 21.34t/m. lEgan & Pricelll996tlEga"n et al.ll2003h . 
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Fig. 2. Examples for the co nstructed light curves in t he V band. 
Data for FF Aql are from iMoffett & BarnesI (Il984t) and from 
iBerdnikovl (l2008h for Y Oph. The resulting light curves are 
drawn with solid lines. 

Our best-fit values are presented in Table |4] and plotted in 
Fig. |3] (black solid curve). This fitted effective temperature is 
only 3 % smaller than the value from iLOSi (6062 + 43 K) at this 
phase of pulsation. The angular di ameter is in excellent agree- 
ment with the 0.86 + 0.17 mas from lGroenewegenl (l2007h . 

We detected a likely IR emission of the order of 2 % (with 
respect to the stellar flux) in PAHl and SiC (see TableO. FF Aql 
presents a more important excess at other wavelengths (8 to 
30/zm). The lower panel in Fig. [3] represents the relative flux 
difference (relative to the model), integrated over the filter band- 
pass. 

We fitted to this star a second component assuming an op- 
tically thin envelope that follows an intensity distribution of the 
form: 



F,x^K,l3B,(Ti) 



(1) 



where B^iTi) is the Planck function at dust temperature T^, is 
the dust opacity and y6 is a parameter directly proportional to the 
dust mass ( Li200l) : 



3 Ma 



^ = 2.1x10--^ 



Fig. 3. Synthetic spectra of the classical Cepheid FF Aql (solid 
line) with the photometric measurements taken from the liter- 
ature. The blue triangles are the points used to fit the SED. 
Our measurements are presented with the red circles while the 
black squares are the other photometric data. On the lower panel 
is plotted the excess flux density relatively to the photospheric 
emission. The green curve represents the flux density model. 



with D the distance of the star in pc and Ma in M©. 

The IR excess seems to rise quite sharply around 10 //m and 
20jum that could be linke d to silicate/carbon features. W e so 
chos e the dust opacity from lOssenkopf & Henning|(ll994i here- 
after |o3. for carbon and silicate grains for a MRN size distri- 
bution ( iMathis et al.|[T977l) . 

We plotted in Fig.|3](the green solid curve) our best flux den- 
sity model with a temperature - 539± 52Kandy6 - 3.0+0. 7x 
10 '''kgm^^. Using the distance from [Benedict et aH (l2007l 
Table [1), this leads to a dust mass ofMd = 1.8 + 0.4 X 10"" Mo 
and to a total mass (gas + dust) of 1.8 + 0.4 x 10"^ Mq (using a 
gas to dust ratio of ~ 100 typical of circumstellar dust). We have 
to mention that some of the additional broadband photometric 
values could be ov erestimated due to th e surrounding interstel- 
lar cirrus emission (iBarmbv et al.ll20in) . 

3.3.2. AX Cir 

As AX Cir has almost the same spectral type than FF Aql we 
considered a surface gravity of log g = 2.0 (note that a change 
of +0.5 only change our following estimated values by 0.7 %). 
Fig, m s hows the SED fit for this Cepheid. B and V photometry 
are fromlB08l (f r om lig ht curves) with a color excess index from 
iTammann et al.l (l2003h . Data were also retrieved from the Deep 
Near Infrared Survey of the Southern Sky (DENIS) for the J and 

bands. Beyond 9 fim data come from MSX, IRC and IRAS. 

Our fitted parameters are presented in Table |4] The di- 
ameter is in agreement with the mean value predicted by 
iMoskalik & Gorynya (2005, hereafter M05, at a 7 % level, con- 
sistent with the angular diameter variation). 

We did not detect IR emission at 8.6 /im and the value 
is consistent with the 9fim measurement from IRC (see 
Table [3). However, it seems to have an infrared excess at 
longer wavelengths. We notice that the IR excess seems to rise 
around 13yL(m and could be a feature of alumina oxide emis- 
sion. Using optical constant for amorphous compact alumina 
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Table 4. Best-fit values for our Cepheid sample. 



Stars 





logs 


6 


LD 








-^cse 


T 

X' 








(mas) 


(K) 


(Ma) 




(K) 




FF Aql 


0.62 


2.05 


0.86 


±0.03 


5890 ± 235 


1.8 ± 0.5 X 10 




539 ± 53 


0.41 


AX Cir 


0.27 


2.00 


0.76 


±0.03 


5911 ± 184 


7.4 ± 5.9 X 10" 




712 ± 61 


0.42 


XSgr 


0.72 


2.00 


1.30 


±0.04 


5738 ±314 


3.0 ± 0.6 X 10" 




703 ± 52 


0.74 


rj Aql 


0.47 


1.80 


1.86 


±0.12 


5431 ±498 


1.0 ± 0.7 X 10" 




545 ± 60 


1.01 


WSgr 


0.48 


1.70 


1.14 


±0.06 


5632 ± 162 


1.7 ± 0.4 X 10- 


-9 


853 ± 55 


1.45 


YOph 


0.73 


1.80 


1.24 


±0.05 


5870 ± 387 


3.9 ± 0.8 X 10" 


-9 


1419 ± 148 


1.49 


UCar 


0.49 


1.20 


0.90 


±0.02 


4823 ± 52 


4.4 ±2.1 X 10" 


-9 


746 ± 94 


0.70 


SV Vul 


0.04 


1.40 


0.76 


±0.01 


5744 ± 144 


3.9 ± 7.4 X 10" 


-8 


620 ± 50 


1.13 


RSct 


0.48 


0.00 


1.74 


±0.06 


4605 ± 119 






1486 ± 335 


1.79 


















772 ± 82 




AC Her 


0.14 


0.50 


0.30 


±0.02 


5711 ±293 






286 ± 32 


0.91 


K Pav 


0.90 


1.20 


1.09 


±0.05 


6237 ± 119 


9.9 ±2.5 X 10- 


10 


695 ± 36 


0.30 



Notes. (f> is the pulsation phase. (?ld and Tetf are the photospheric fitted parameters while 7^^ is the CSE fitted parameter. M^se is the total mass 
(dust + gas) of the CSE estimated using a gas to dust ratio of ~ 100. For R Set, the first line denotes the compact component. The^- has been 
estimated from the fit of the photospheric flux. See the text for the references of the effective gravity. 



dBegemann et aljri997h . we computed from a simple model 
of h omogeneous spheres in the Rayleigh limit of small parti- 
cles ("Bohren & Huffmanl \mf), with a size distribution from 
Mathis et al., (19771 with particle sizes from 5 nm to 250 nm as 
094I) and a density p = 2.5 g cm . We then used this dust opac- 
ity in the model of Eq.[T|to fit a second component (green curve 
in Fig. |4]i, assuming that the envelope is mainly composed of 
amorphous alumina. We obtained /3 - 1.6 + 0.4 x lO"'** kg m"-. 
Table |4] shows the fitted temperature and the total mass of the 
CSE estimated using a gas to dust ratio of ~ 100. However, a 
larger photometric dataset is necessary to better constrain the 
CSE parameters and its dust composition. 




2. 5. 

?teff (fim) 

Fig. 4. The same as Fig. |3] but for AX Cir. 



3.3.3. X Sgr 

We retrieve photometric data from the B, V dereddened magni- 
tudes of Kervella et al. (,2004a hereafter iK04h and from lM84l 
(from light curves at an interm ediate phase (h - 0. 72). The other 
photometric data come from iFeast et al ] (i2008l from J, H, K 



light curves, hereafter|F08|), IRC (9 //m), IRAS (12|/m) and MSX 
(8.28, 12.13 and 14.65 jum). These data are plotted inFig.|5] 

For the same surface gravity previously used {\ogg - 2.0) 
we found an angular diameter (see Table |4|i th at is 13% and 
3cr sma ller than the mean diameter measured bv lKervella et al.l 
(l2004bi 1.47 + 0.04 mas, see Table [B. How ever, the ampHtude 
of the pulsation is ~ 9 % in diameter (lM05h . Conversely our es- 
timate is in a greement with 0LD - 1-31 +0.1 2 mas assessed from 
the parallax (Benedict e t al.ll2007h a nd the linear diameter at this 
phase of pulsation (FOSj). However lK04l used a limb darkened 
diameter to model their data and the presence of an extended 
emission could overestimate the angular diameter. 

We detected a likely IR emission of the order of 5-15 % 
(with respect to the stellar flux, see Table |3]l. We fitted the same 
flux density model previously used, with the dust opacity from 
094 as it seems to have some silicate/carbon features. We ob- 
tained p - 5.6 + 0.9 X lO^'^kgm"^. The other parameters are 
presented in Table |4] 




2. 5. 
^Leff (|im) 

Fig. 5. The same as Fig. [3] but for X Sgr 



10.0 20. 
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3.3.4. J] Aql 

Optic al and nea r-infrared pho tometry from 



K04l {B, V, J, H, K), 



IM84I {B, V) and iBarnes et al.l {J, H, K, Il997h were used for the 
SED (from light curves at an intermediate phase = 0.47). 
For the longer wavelengths we retri eve the IRAS fluxe s and the 
Spitzer/IRS spectra (from 5 to 35 ^m. lArdila et al.ll2010l) . The re- 
sult is plotted on Fig. |6] The IRS spectra is plotted in purple. 
This spectra was acquired at a phase - 0.45 very close to our 
intermediate phase of our observations and therefore we did not 
add uncertainty due to the phase mismatch. In the lower panel of 
Fig. |6]we just plotted the spectra for some valu e (bin of IQfim) 
for clarity. The Spitzer flux measurements from lMlOl were also 
added, but are significantly lower than the stellar SED derived 
from visible and near-IR. A possible explanation of this incon- 
sistency could be due to the different observing mode used for 
J] Aql with respect to the other Cepheids. tj Aql was observed 
in full-frame mode while the other stars used in this paper were 
observed in subarray mode. In the full frame mode the images 
of this bright star are heavily saturated. The PSF-fitting method 
used in |M13 could be incorrect, leading to an underestimation of 
the photometry. For this reason, we decided not to include them 
in the fit. 

We found 0ld (see Table IHl in good agree ment with the 
measured diameter from iKervella et all ( l2004bt 0ld — 1.87 + 
0.03 mas) at this phase and als o with the mean angular diam- 
eter from lGroenewegenl(l2008l p - 1.76 + .09 mas). At an 
intermediate phase (<p ~ 0.47) iLuck & Andriev skv (2004!) give 
a surface gravity of \ogg ~ 1.8 and an effective temperature of 
about Teff = 5508 + 40 K that is only 1 % larger than our esti- 
mated value. 

We can see from the IRS spectra that it is likely that this 
Cepheid has a small IR emission, starting with a relative excess 
(still with respect to the photosphere) of 4.7 + 1 .5 % at 5.3 yum 
to 9.2 + 1.8 % at 34.7 pm. We also detect in our VISIR filters an 
excess with the same order of magnitude (see Table[3]l. 

No strong features of silicate seem to be present in the spec- 
tra. We so made the hypothesis that the envelope is mainly com- 
posed of amo rphous carbon. Using optical constant for amor- 
phous carbon dPreibisch et al.l[T993h . we also computed k,] from 
a simple model of homogeneous spheres in t he Rayleigh hmit of 
small particles, with a size distribution from lMathis et al.l (1 19771 
with particle sizes from 5nm to 250 nm as I094I) and a den- 
sity p - 2.0 gcm"^. We then used this dust opacity in Eq. [T] 
to fit the infrared excess (green curve in Fig. |6]l. We obtained 
j3 - 1.6 + 0.4 X 10"'^ kg m"- that gives the total dust mass pre- 
sented in Table|4]with the dust temperature. 

3.3.5. WSgr 

The B and V photometry are from lK04l and lB08l (at an interme- 
diate phase (p - 0.48). The other irradiances ar e from DENIS 
iJ,Ks), Spitzer (3.6, 4.5, 5.8, 8 and 24 ^um from lMlOl) . IRC (9 
and l8pm),MSX (8.28, 12.13, 14.65, 21.34//m) and the 12pm 
photometry from IRAS. The SED is plotted in Fig.|2] 

Our fitted effective temperature (presented in Table HJi is 
in agreement at a 2 % lev el with T^ff - 5535 + 51K from 
iLuck & Andrievskv! (|2004|) around the same phase with an ef- 
fective gravity from the same author of 1.7. We found an an- 
gular diameter that is about 15 % and 2cr smaller than the one 
measured by interferometry bv |K04| (1.31 ± 0.04 mas) at this 
phase of pulsation. Conversel y our estima t e is in agreement 
with the mean diameter from iBersier et al ] (Il997h who found 
&LD - 1.17 + 0.11 mas based on photometry. However lK04l used 




0.5 1 .0 



Xeff (nm) 



Fig. 6. The same as Fig.|3] but for rj Aql. The measured spectra 
were superimposed in purple. 



a limb darkened diameter to model their data and the presence of 
an extended emission could overestimate the angular diameter 

Spitzer' s valu es are consistent with the blackbody radiation 
adjusted bv iMlOl However our VISIR photometry shows an ex- 
cess of ± 15 % (see Table [3]). The same trend is observed from 
IRC (at 9pm) and MSX (at 8.28, 12.13 and 14.65 /zm). This 
could be an evidence of a particular dust composition. As no 
IR emission is detected around 20pm, we rejected a silicate 
dust composition and we so assumed an envelope mainly com- 
posed of alumina oxide. We used the same dust opacity model 
as AX Cir in Eq. [T] The resulting curve is plotted in Fig. [T] and 
the fitted parameters are listed in Table |4] 




Fig. 7. 



>^eff (lirn) 

The same as Fig. [3] but for W Sgr 



3.3.6. YOph 

We retrieved th e B, V data fromlBOSi a nd lM84 The J, H, K pho- 
tometry is from lLanev & Stobld7T992h (al so from light curves at 
an intermediate phase (p - 0.73) and K04l The plot of the spec- 
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tral energy distribution is shown in Fig. |8] T he oth er fluxes are 
from Spitzer (3.6, 4.5, 5.8, 8 and 24 ^um from lMlOh . IRC (9 yum) 
and the llfim and 25 fim photometry from IRAS. 

The best fitted parameters are presented in Table |4] for a 
found logg ~ 1. 8 dLOSh . |L08| also give an effective tempera- 
ture at this intermediate phase of T^ff = 5800 + 148 K that is 
in good agreement with our value. Our fitted angular diameter 
and the value 0ld - 1-24 ± 0.01 mas at <p ~ 0.73 measured by 
iMerand et all ( l2()07h is also good. 

We detected a likely infrared emission from our photomet- 
ric measurements (s ee Table O. This result is consistent with 
iMerand et al.l (l2007h where a CSE has been found around this 
star in the K band with a relative contribution of 5.0 + 2.0 %. 
As some features of silicate/carbon seems to be present around 
lOfim and 20 /im, the model of Eq. [1] was fitted with the dust 
opacity from 094. The best-fitted parameters are presented in 
Table H] Y Oph seems to have a hotter circumstellar envelope 
that could be located close to the star and heated by its radia- 
tions. 




2. 5. 

Fig. 8. The same as Fig. |3] but for Y Oph. 



3.3.7. U Car 



We selected the B, V photometry from lB08l and 
ICoulson & Caldwelll (Il985h (from light curves at an mter- 
mediate phase = 0.6 3). We also used the /, H, K values from 
iLanev & StoblS (119921 hereafter L92, also from light curves) 
and mid-IR photometric measurements from IRC (9 and 1 8 //m), 
from MSX (8.2 8, 12. 13, 14.65 yum) and Spitzer (3.6, 4.5, 5.8, 8 
and 24 /im from lMlOh . In Fig.|2]is plotted the SED. 

Our fitted angular diam eter and an effectiv e temperature with 
a given log g ~ 1.2 from iRomaniello et al.l (I2008D (given at a 
phase of 0.49 but as said previously the broadband photometry 
is almost insensitive to the effective gravity changes) are listed 
in Table |4] Our estimate of the diameter is in agreement a t a 5 % 
level with the mean angular diameter from lGroenewegenI (l2008l 
0.94 + 0.05 mas). However our estimated effective temperature 
is not consistent with the me an value Tf.g = 5980 K e stimated 
from the iron abundances by Romaniello et al.' (200S). From a 
color-temperature relation (Fry & Carney 1999), the tempera- 
ture of the star should be T^g ~ 50 00 K making suspicious the 
value from lRomaniello et al.l (l2008h . We did not find other mea- 



sured temperatures in the literature to verify our estimate but we 
are confident in our value since all photometric data in the SED 
(< 8 jj.m) look consistent with this effective temperature. 

We have detected with VISIR a significant IR emission (see 
Table |3]l. This is also visible at larger wavelength with an in- 
frared excess of about 20-30 %. We also fitted the flux density 
model (Eq. [U with the dust opacity from '094l The results are 
listed in Table |4] We should notice that these values could be 
affected by the i nterstellar cirrus bac kground emission in the re- 
gion of this star (iBarmbv et al.ll20ll1) . 




0.5 1.0 2. 5. 

Fig. 9. The same as Fig. [3] but for U Car 



3.3.8. SV Vul 



ICO 20. 



The B, V photometry were retrieved from lBOS] and lM84l while 



the J, H, K magnitud e s from light curves of lBarnes et al. 
and iLanev & Stobld (Il992l) . Other infrared data are from IRC 
(9 and 18/im),M5X(8.28, 12.13, 14.65 yum) and (12jum). 
The spectral energy distributi on is presented i n Fig.fTOl 

We took logg ~ 1.4 from lKovtvukh et al.1 (120051) . These au- 
thors give an effective temperature around the same phase of 
TeSf = 5977 + 32 K. We found a temperature that is only 4 % 
and 1.4cr smaller The an gular diameter we fou nd is 5 % smaller 
than the mean value from lGroenewegenI ( l2008i 0.80 + 0.05 mas), 
but within Icr. 

We have detected an excess in PAHl and in SiC filters (see 
Table O. This excess is also visible at other wavelength up to 
26.7 + 6.6 % at 18//m with IRC. This indicates the presence of 
circumstellar material. 

We fitted the same model previously used (Eq. [U for the 
wavelengths larger than 3 pm and still with the dust opacity from 
094, The fitted parameters are listed in Table H) The uncertainty 
in the total mass estimate of the envelope is due to the uncer- 
tainty on the parallax measurement (Table[T]). More photometric 
data points are necessary to better constraint this SED, measure- 
ments at least between 3 and 8 pm should be useful. 

3.4. Spectral energy distribution of type II Cepheids 

In this section we use the same methods and models for type II 
Cepheids than for classical Cepheids. 
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0.5 1.0 2. 5. 

Fig. 10. The same as Fig. [51 but for SV Vul. 
3.4.1. R Set 

This very irregular pulsator with a large amplitude has been ex- 
tensively observed and is known to have an infrared exc ess (see 
e.g.iG oldsmith et al. 1987; Giridhar et al. 2000; De Ruvt er et al.l 
l2005h . Its irregular pulsation and its large amplitude make it 
a special case as R V Tauri stars. The B, V,J,H,K photometry 
were retrieved from IShenton et al.l (1 19941) (at phase 0.48). We 
also added values from IRC (9 and 18jum), IRAS (12yum) and 
from MSX (8. 28, 12.13, 14.65 21.34 urn) E{B-V) and B-V w ere 
obtained from lTaranova et alJ (1201 Ol) and lMvers et alJ (1200 ih re- 
spectively. The SED is plotted in Fig.fTTl 

We fitted the stellar component (i.e. A < 3 nm) wit h an ef - 
fective gravity of log^ = 0.0 given bv iGiridhar et all (l2000l) . 
As expected we also detected a strong IR emission at our wave- 
lengths (see Table [3]l. The fitted photospheric curve is shown 
in Fig. [TT] (solid black curve) and its paramete rs in Table HI 
Our te mperature is consistent at a 2 % level with IGiridhar et all 
(I2OOOI) wh o gave T^ff = 4500 K (at (p = 0.44 using our 
ephemeris). IShenton et al.l (1 19941) estimated a diameter of about 
1 .79 ± 0.59 mas (taking d - 43 1 pc from Hipparcos) and our 
value is only 3 % smaller 

For the circumstellar component, we assumed that the enve- 
lope is optically thick and the observed flux density can be fitted 
by the following equation: 



F,i^AB,{Ta) 



(2) 



where A is related to the solid angle and the emissivity of the 
envelope. 

By looking at the SED our values seem to be lower than the 
other measurements and it seems to have two different intensity 
distributions. As the effective resolution of VISIR is larger than 
that of the AKARI and IRAS telescopes in the band, VISIR only 
detected a compact component located close to the star. We so 
decided to first fit only our values as a compact component and 
then to fit a larger component using only the IR data from AKARI 
and IRAS. The result are presented in Fig.[TT] The solid red curve 
is the possible warmer compact component for which our best 
fitted values were T = 1486 + 335 K and A = 4.65 + 0.83 mas^. 
For the larger component (the green curve) we found T - 772 + 
83 K and A = 1 1 .2 + 1 .25 mas^ (Tableg]). This latter temperature 
is in the sa me range than the dust temperature T = 800 + 50 K 
modelled bv lTaranova et al.l (I2OIOI) . 



Fig. 11. Synthetic spectra of the type II Cepheid R Set (solid 
line) with the photometric measurements taken from the litera- 
ture. In blue are the points used to fit the SED. Our measure- 
ments are presented in red. On the lower panel is plotted the 
excess flux density relatively to the photospheric emission. The 
red curve represents the flux density model for a compact com- 
ponent while the green curve denotes larger one. 



3.4.2. AC Her 

AC Her has been studied extensively for its strong IR excess 
due to thermal emission from circumstellar dust grains (see 
e.g. iDe Ruvter e"tal1l2005l). We ret r ieved fi, V, J, H, K band pho- 
tometric flux from Ishenton et all ( 1 19921) (at phase 0.14). The 
longer wavelengths are data from IRC (9 and 18//m), IRAS 
(12, 25, 60, lOO^um ). E (B-V) and B-V we re retrieve from 
iTaranova et"ai] (|20IO|) and iMvers et alj (1200 Ih respectively. We 

plotted the SED in Fig.lT2l 

As in iDe Ruvter et al.l (l2005h an infrared excess is clearly 
visible. Our values lie on the maximum of the spectral energy 
distribution and also show a strong contribution from the cir- 
cumstellar environment (Table [3]). We fitted a two component 
spectrum corresponding to the photosphere (still using Kurucz 
model) and the circumstellar emission. The best fitted parame- 
ters found for the photospheric emission (black soUd curv e in 
Fig. [nil with logg = 0.5 from IVan Winckel et alJ (Il998l) are 
in Table |4] The temperatu r e is in agreement with the one es- 
timated bv ITaranova et al.l (l20ICt mean value: T^ff ~ 5400 K 



with Tq - 5 8 00 K) . The angular diameter is consistent with 
IShenton et aD (1 19921 0ld = 0.31 mas with R = 23.6 /?o and 
d - 750 pc). The second component was fitted with the model 
of Eq.|2l assuming an optically thick medium (solid green curve 
in Fig.lTSI). The fitted temperature is listed in Table |4| with A = 
92.0 ±6.9 mas~. Our esti mated temperature is larger than the one 
from lGielen et al.l (l2007h who found a blackbody temperature of 
170 K using a disc model, his model also gave Ri^ = 50 mas and 
Rout - 428 mas for the disc (with d = 1400 pc). However the 
author assumed that the morphology of circumstellar material 
is a disc. Other values are found from a dust shell model as- 
sumed by iDe Ruvter et al. ( 2005h who found R-^ - 22 mas and 
Rout - 397 mas. lAlcolea & Buiarraball (Il99ll) also fitted the ex- 
tended emission with a shell model and found R^n - 76 mas and 
Rout - 3330 mas with a dust temperature of T - 360 K. We can 
also notice that AC Her was also studied using high-resolution 
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imaging with adaptive optics bv lClose et alJ(l2003h who did not 
detect any extended mid-infrared structure greater than 200 mas. 




0.5 1.0 2. 



5. 10.0 20. 



50. 100.0 



Fig. 12. The same as Fig. (TT] but for AC Her 
3.4.3. K Pav 

We retrieved the V, J, H, Ks photometric light curves from iFOSl . 
We used E{B - V) end B-V from F08. The mid- and far-infrared 
data were obtained from IRC and IRAS. We plotted the SED in 

Fig.iia 

W e chose a mea n effective gravity of log^ - 1.2 from 
iLuck & Bondl (Il989l) . Our fitted parameters are presente d in 
Table |4] Our diameter estimate is in good agreement with |F08| 
(1.04 + 0.04 mas (p - 0.90, with^f = 204 + 6 pc). However we 
found an eff ective temperature 9 % larger than the mean value 
estimated bv iLuck & Bondl ([T9891 T^ff ~ 5750 K at = 0.94) 
based on metal abundance analyses. It has to be noted that this 
star is classified as peculiar type 11 Cepheid with distinctive light 
curve and brig hter than nor mal type 11 Cepheid with the same 
period (see e.g. lMatsunaga e t al. 2009). 

We detected with VISIR an IR excess of the order of 20%. 
This excess tends to increase with wavelength leading to the hy- 
pothesis of an extended circumstellar envelope. We fitted to the 
IR photometric data the model of Eq. [T] assuming the envelope 
is optically thin. As some features are present around lO^um and 
20 /im, we also assumed a co mpos ition of silicate/carbon grains. 
We used the dust opacity of l094l to derive the total dust mass, 
assuming a gas to dust ratio of 100. The best fitted parameters 
are listed in Table H) 

3.5. Spatially resolved emission 

We search for spatially extended emission using a Fourier 
technique, similar in its principle to the calibration tech- 
nique used in long baseline interfe rometry. This method 
was already used and validat e d by iKervella et alj (|2009|) : 
iKervella & Domiciano de Souzal (120061) . The principle is to di- 
vide the Fourier transform modulus of the image of the Cepheid 
(^cep) by that of the calibrator image (/cai): 



*P(v.,v,) = 



2. 5. 10.0 20. 

Xeff (nm) 

Fig. 13. The same as Fig. (TT] but for k Pav. 



where the hat symbol denote the Fourier transform, {x,y) the 
sky coordinates and (Vj, v,.) the angular spatial frequencies. This 
equation is related to interferometric observations that provide 
measurements of the Fourier transform of the intensity distribu- 
tion of the observed object (Van Cittert-Zernike theorem). 

We then compute the ring median of i.e. the median for 
a given spatial frequency radius v over all azimuths (where = 
Vj. + v^). The function ^(v) obtained is equivalent to a visibility 
in interferometry. The error bars on ^' were estimated by the 
quadratic sum of the dispersion of the PSF calibrators Fourier 
modulus over the night and the rms dispersion of the calibrated 
function over all azimuth directions for each spatial frequency. 
A deviation from a central symmetry will not be detected and 
any departure will be included in the error bars. 

Defining a model of a point-like star surrounded by a 
Gaussian shaped CSE and taking its Fourier transform, it 
is possible to retrieve the CSE i ntensity distribution. This 
type of model was already used by Kerv ella et al.l (l2009h and 
iKervella & Dom iciano de Souza. (.2006.) with the *P(y) function 



>P(v,p,,a,,) = 



/* + /cs 



1 + ai 



I + exp 



41n2 



where the Gaussian CSE is defined with a FWHM and a rel- 
ative flux - fcsei'l)/ f*{^), i-e the ratio of the flux of the en- 
velope to the photospheric flux. We set = 1 since the star is 
unresolved by VISIR. 

We applied the fit using a classical minimization to all the 
final images. We did not detect spatially resolved emission for 
FF Aql, I] Aql, U Car (in the SiC filter), SV Vul, R Set, AC Her 
and K Pav. For the other Cepheids that show a resolved compo- 
nent, the fitted parameters are presented in Table |5] and the ^' 
function are plotted in Fig.fTSl 

For the stars with an undetected envelope, an upper limit of 
~ 265 mas can be set for the extension of the CSE based on the 
telescope resolution. 

The fact that we detect an extended emission around AX Cir 
and not in the previous section (Sect. 13.3b can be explain by the 
lack of J,H,K light curves to fully constrain the SED at our 
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Fig. 14. Sample of the final VISIR images for our Cepheids and their calibrator The scale is logarithmic for all images. 



we derive a minimum mass loss rate M ~ 7 x 10 M^yr ' 
for AX Cir and W Sgr, and M ~ 2 x IQ-^^ M^yr^ for X Sgr, 
Y Oph and U Car Thes e values are comparable with the pre- 
dicted mass loss rate from lNeilson & Lested (l2008l ranging from 
10~'" t o lO^ ^Mpyr^') and the measured mass loss rate from 
iDeasvl (119881 ranging from 10"'° to IQ-^ M^yr-^). 

The Fourier technique gives IR excesses that are very compa- 
rable to those derived from the SED fitting presented in Sect. 13.31 
(except for AX Cir). This convergence of two independent anal- 
ysis methods gives us good confidence that the derived envelope 
parameters are reliable. 

4. Discussion 

Classical and type II Cepheids are two distinct classes of pul- 
sating stars. The formers are know to be intermediate mass stars 
with regular pulsation periods, evolving as post-main sequence 
stars in a core He-burning phase. Type II Cepheids are gener- 
ally associated to lower mass stars with irregular pulsations and 
undergoing the post-core He-burning evolution. Therefore the 
CSEs surrounding these two classes correspond to different evo- 
lutionary stages. 

R Set and AC Her have been observed extensively and are 
known to have strong IR excesses linked to the extended CSEs 
which are interpreted as being a relic of their strong dusty mass 



Table 5. Fitted parameters of the 'V{v,px, a\) function. 



Stars 


Filter 


P 




a 






(" 


) 


(%) 


AX Cir 


PAHl 


0.69 ± 


0.24 


13.8 ±2.5 


XSgr 


PAHl 


0.99 ± 


0.23 


7.9 ± 1.4 




PAH2 


0.99 ± 


0.44 


15.2 ±3.7 




SiC 


0.81 ± 


0.53 


8.9 ±3.5 


WSgr 


PAHl 


1.14 ± 


0.39 


3.8 ±0.6 




PAH2 


1.19± 


0.37 


9.1 ± 1.5 




SiC 


1.03 ± 


0.50 


8.3 ± 2.4 


YOph 


PAHl 


0.71 ± 


0.12 


15.1 ± 1.4 




PAH2 


1.02 ± 


0.52 


7.5 ±2.3 




SiC 


0.54 ± 


0.46 


6.2 ±4.1 


UCar 


PAHl 


0.74 ± 


0.10 


16.3 ± 1.4 



phase of pulsation. The uncertainties in these bands (0.14 mag 
in K) are large enough to prevent a detection of a IR excess. 

From the distance p estimated with this technique, it is possi- 
ble to assess a lower limit for the mass loss rate by computing the 
time t required by a stellar wind to re ach this distance. This c om- 
putation has already been applied bv lMarengo et al.l (1201 Obi) for 
the Cepheid 6 Cep. Us ing the escape velocity Vesc ~ 100 km s ' 
(IWelch & Duricl[l988h as the minimum wind speed, we found 
/ ~ 12-23 yr With the total dust mass estimated in Sect. 13.31 
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U Car 




PAH1 

0.0 0.5 1.0 1.5 2.0 

Spatial frequency (arcsec"''' 

Fig. 15. Fit of a Gaussian CSE + unresolved point source model to the ^l* functions of our Cepheids. The dashed curve is here for a 
reference as an unresolved star. 
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loss on the asymptotic giant branch. Different models have been 
applied to constrain the morphology of this surrounding mate- 
rial. De Ruyter et al. (2005) interpreted the spectral energy dis- 
tribution of these stars using a dust shell model. They found an 
inner radius Ri^ - 12.5 AU and an outer radius i?out = 224 AU 
for AC Her and Rin = 13.4 AU and Rp^t ^ 5200 AU for R Set. 
From a disc model applied to AC Her. lGielen et al.l (l2007h found 
Rin = 35 AU and a smaller external radius Rput = 300 AU. For 
this star these results are in contradiction with lClose et al. I (l200l 
who excluded any extended emission larger than 75 AU from 
high Strehl ratio adaptive optics images in the mid-lR domain. 
We also not detected in our filters with VISIR extended emission 
larger than 7? = 100 AU. 

The case of R Set is peculiar a nd this star is often classified 
as an exception by several authors. lAlcolea & Buiarrabail ( ll991h 
used a two shell model to interpret both the mid-IR and far-IR 
data and they estimated a temperature for the largest grains of 
T - 815 K. The estimated radii of the shells were Rin,\ - 30 AU, 
^out,i = Rin,! = 5 720AU and R^^t.i = 12 000 AU (using d = 
43 1 pc) where indexes 1 and 2 denote the inner and outer shell 
respectively. We possibly detected the emission of the inner shell 
from our VISIR photometry but no spatially extended emission 
larger than = 100 AU has been detected. 

K Pav has been recently classified as a peculiar type II 
Cepheid by some authors because of its distinctive light curve. 
It is significantly brig hter than normal type I I Cepheids with the 
same period (see e.g. lMatsunaga et ani2009h . This star is not at 
the same stage of evolution as R Set and AC Her since it is as- 
cending in the HertzsprungRussell diagram the blue horizontal 
branch to the asymptotic giant branch. During this process the 
star undergoes changes both in its core and envelope that could 
lead to mass loss via pulsation and/or shock mechanisms. An 
active past and/or on-going mass loss may explain the large IR 
excess we have detected. 

With classical Cepheids we are on different scales since 
the envelope that have been detected up to now only have a 
spatial extension of a few stellar radius. Only a small sam- 
ple of classical Cepheids is known to host a circumstellar en- 
velopes ({ Car, Y Op h, RS Pup, Polaris, 6Cep, S Mus, GH Lup 
T M on and X Cvg. iKervella et alJ 120061: iMerand et alJ l2006l 



ence of sho ck waves at certain pulsation phases (iNardetto et al.l 
I2006ll2008h . 



1 M on ana A uvg. | Kerveiia et aiJ izuuot iMerana et aiJ izuua. 
l2007HBarmbv et alJl201 ll) . With this paper we extend this sam- 
ple with FF Aql, tj Aql, W Sgr, U Car and SV Vul as probable 
IR excess Cepheids. The origin of these CSEs is not well un- 
derstood. Their presence could be linked to mass loss from the 
star. It has been suggested that the I R excess is caused b y dust 
formed in a wind from the Cepheids dKervella et alJ2006l) . From 
a rad iative-driven wind model including pulsation and shock ef- 
fects 'Neilionl&Leitei (120081) concluded that radiative driving is 
not sufficient to account for the observed IR excesses and pro- 
posed that the mass loss could be driven by shocks generated in 
the atmosphere by the pulsation of the star. Other observational 
evidences were provided with IR excess detections from IRAS 
observations bv iDeasvl (Il988l) who estimated mass loss rates 
ran ging from 10~ ' ° to 10 '^ Mq yr"'. 

iMerand et al.l (l2007h showed for the classical Cepheids a 
likely correlation between the pulsation period and the CSE flux 
(relative to the photosphere) in the K band. We plotted in Fig. [16] 
the relative CSE flux versus the period at 8.6^(m. It seems that 
at 8.6 A<m the lon g er per iods have larger excesses, as concluded 
by IMerand et al.l (|2007|) in the K band. Assuming this excess 
is linked to a mass loss phenomena, this correlation shows that 
long period Cepheids have a larger mass loss than shorter period, 
less massive stars. Such a behavior could be explained by the 
stronger velocity fields in longer period Cepheids, and the pres- 
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Fig. 16. Measured relative CSE fluxes at 8.6//m around 
Cepheids as a function of the pulsation period. 



5. Conclusion 

We presented new thermal IR photometry and spectral energy 
distributions of 8 classical Cepheids and 3 type II Cepheids. We 
detect extended emission around 10 stars in the thermal infrared 
domain, of which 3 were previously known. The sample of clas- 
sical Cepheids with CSEs is extended with five additional stars 
to a total of 9 stars. This confirms that the presence of circum- 
stellar material around classical Cepheids is a widespread phe- 
nomenon. 

A correlation is also found between the pulsation period and 
the CSE flux in the thermal IR and confirms the same correla- 
tion found in the K band by IMerand et"al] ( l2007h . It thus con- 
firms that longer period Cepheids have apparently larger mass 
loss rates than shorter periods. Apart from their probable impor- 
tance in the evolution of Cepheids, the existence of these CSEs 
may also impact Cepheid distance measurements in the infrared 
domain, in particular with the JWST. 

We also point out the need for more optical and near-infrared 
Cepheid light curves. This is the main source of uncertainty of 
the spectral energy distributions at a given phase of pulsation. 
A dedicated J, H, K photometric survey of few month would be 
sufficient for a good coverage in phase and would favour the 
mid- and far-infrared excess detection. 
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Table 2. Log of our VISIR BURST mode observations. 



MJD 





Stai- 


Filter 


DIT 

(ms) 




seeing 
(") 


AM 


# 


54 610.035 




HD 89682 


PAHl 


16 


22 500 


1.2 


1.26 


1 


54 610.042 




HD 89682 


PAH2 


8 


48 000 


1.1 


1.27 


2 


54 610.056 


0.62 


UCar 


PAHl 


16 


22 500 


1.5 


1.30 


3 


54 610.064 


0.62 


UCar 


PAH2 


8 


48 000 


1.1 


1.32 


4 


54 610.081 




HD 98118 


PAHl 


16 


22 500 


1.1 


1.31 


5 


54 610.088 




HD 98118 


PAH2 
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48 000 


1.1 


1.35 


6 


54 610.104 


0.65 


XSgr 


PAHl 


16 


22 500 


1.0. 


1.61 


7 


54610.111 


0.65 
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PAH2 
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24 000 


1.0 
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8 


54 610.126 


0.71 


YOph 


PAHl 


16 


22 500 


1.2 


1.63 


9 
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0.71 
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PAH2 


8 


48 000 


1.8 
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10 
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PAHl 
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11 
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48 000 
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48 000 
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48 000 
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T] Aql 
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48 000 
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0.02 


SV Vul 
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HD 203504 


PAH2 


8 


48 000 


0.9 


1.63 


26 


54 610.349 


0.14 


AC Her 


PAHl 


16 


22 500 


0.8 


1.58 


27 


54 610.357 


0.14 


AC-Her 


PAH2 


8 


48 000 


0.9 


1.60 


28 


54 610.373 


0.71 


YOph 


PAHl 


16 


22 500 


0.8 


1.32 


29 


54 610.380 


0.71 


YOph 


PAH2 


8 


48 000 


0.7 


1.36 


30 


54 610.395 


0.42 


rj Aql 


PAHl 


16 


22 500 


0.8 


1.15 


31 


54 610.402 


0.42 


rj Aql 


PAH2 


8 


48 000 


0.9 


1.17 


32 


54 610.417 




HD 196321 


PAHl 


16 


22 500 


0.8 


1.11 


33 


54 610.424 




HD 196321 


PAH2 


8 


48 000 


0.7 


1.12 


34 


54611.016 




HD 89682 


PAHl 


16 


22 500 


0.7 


1.22 


35 


54 611.024 




HD 89682 


SiC 


20 


18 000 


0.8 


1.24 


36 


54 611.035 


0.64 


UCar 


PAHl 


16 


22 500 


0.8 


1.27 


37 


54 611.042 


0.64 


UCar 


SiC 


20 


18 000 


0.8 


1.28 


38 


54 611.054 


0.26 


AXCir 


PAHl 


16 


22 500 


0.8 


1.40 


39 


54 611.061 


0.26 


AXCir 


SiC 


20 


18 000 


0.8 


1.38 


40 


54 611.073 




HD 124294 


PAHl 


16 


22 500 


0.7 


1.06 


41 


54611.081 




HD 124294 


SiC 


20 


18 000 


0.7 


1.05 


42 


54 611.093 


0.27 


AXCir 


PAHl 


16 


22 500 


0.8 


1.32 


43 


54611.101 


0.27 


AXCir 


SiC 


20 


18 000 


0.7 


1.31 


44 


54611.112 


0.79 


XSgr 


PAHl 


16 


22 500 


0.8 


1.50 


45 


54611.119 


0.79 


XSgr 


SiC 


20 


18 000 


1.1 


1.44 


46 


54611.131 


0.54 


WSgr 


PAHl 


16 


22 500 


0.8 


1.43 


47 


54611.139 


0.54 


WSgr 


SiC 


20 


18 000 


0.8 


1.38 


48 


54611.151 




HD 161096 


PAHl 


16 


22 500 


0.8 


1.49 


49 


54611.158 




HD 161096 


SiC 


20 


18 000 


0.8 


1.43 


50 


54611.170 


0.76 


YOph 


PAHl 


16 


22 500 


0.7 


1.27 


51 


54611.177 


0.76 


YOph 


SiC 


20 


18 000 


0.7 


1.23 


52 


54611.189 


0.48 


RSct 


PAHl 


16 


22 500 


0.6 


1.39 


53 


54611.196 


0.48 


RSct 


SiC 


20 


18 000 


0.6 


1.34 


54 


54 611.207 




HD 161096 


PAHl 


16 


22 500 


0.6 


1.21 


55 


54611.215 




HD 161096 


SiC 


20 


18 000 


0.6 


1.19 


56 


54 611.227 


0.90 


K Pav 


PAHl 


16 


22 500 


0.7 


1.45 


57 


54 611.235 


0.90 


K Pav 


SiC 


20 


18 000 


0.7 


1.43 


58 


54 611.248 


0.53 


rj Aql 


PAHl 


20 


18 000 


0.7 


1.38 


59 


54 611.255 


0.53 


77 Aql 


SiC 


25 


14 400 


0.7 


1.33 


60 
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Table 2. continued. 



MJD 





Star 


Filter 


DIT 




seeing 


AM 


# 










(ms) 




(") 






54 611.273 




HD 133774 


PAHl 


16 


22 500 


0.8 


1.35 


61 


54 611.280 




HD 133774 


SiC 


20 


18 000 


0.8 


1.40 


62 


54611.310 




HD 161096 


PAHl 


16 


22 500 


0.7 


1.21 


63 


54611.317 




HD 161096 


SiC 


20 


18 000 


0.7 


1.23 


64 


54 611.340 


0.15 


AC Her 


PAHl 


16 


22 500 


0.8 


1.54 


66 


54 611.347 


0.15 


AC Her 


SiC 


20 


18 000 


0.8 


1.57 


67 


54 611.365 


0.04 


SV Vul 


PAHl 


16 


22 500 


0.7 


1.64 


68 


54 61 1.3 /2 


A A/1 

0.04 


h\ Vul 


iiiC 


OA 

20 


1 O AAA 

is 000 


A O 
0.0 


1 c 


69 


54 611.384 




HD 203504 


PAHl 


16 


22 500 


0.7 


1.42 


70 


54 611.392 




HD 203504 


SiC 


20 


18 000 


0.8 


1.41 


71 


54 611.403 


0.62 


FF Aql 


PAHl 


16 


22 500 


0.6 


1.66 


72 


54611.411 


0.62 


FF Aql 


SiC 


20 


18 000 


0.7 


1.73 


73 


54 611.422 


0.56 


rj Aql 


PAHl 


16 


22 500 


0.7 


1.24 


74 


54 611.430 


0.56 


rj Aql 


SiC 


20 


18 000 


0.7 


1.27 


75 



Notes. MJD is the Modified Julian Date at the start of the exposures of the target, (j) is the phase of the Cepheid. DIT denotes the Detector 
Integration Time for one short exposure image. N represents the total number of frames. The seeing is measured in the visible by the DIMM 
station. AM denotes the airmass. 



